Hydromagnetic Instability in Differentially Rotating Flows 
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We study the stability of a compressible differentially rotating flows in the presence of the magnetic 
fleld, and we show that the compressibility profoundly alters the previous results for a magnetized 
incompressible flow. The necessary condition of newly found instability can be easily satisfied in 
various flows in laboratory and astrophysical conditions and reads BsB^Q! 7^ where Bs and 
B^p are the radial and azimuthal components of the magnetic field, Q.' = dQ/ds with s being the 
cylindrical radius. Contrary to the well-known magnetorotational instability that occurs only if 
n decreases with s, the instability considered in this paper may occur at any sign of Jl'. The 
instability can operate even in a very strong magnetic field which entirely suppresses the standard 
magnetorotational instability. The growth time of instability can be as short as few rotation periods. 

PACS numbers: PACS numbers: 47.20.-k, 47.65.+a, 95.30.Qd 



INTRODUCTION 

Instabilities caused by differential rotation of a mag- 
netized gas may play an important role in enhancing 
transport processes in various astrophysical bodies and 
laboratory experiments. It is well known since the clas- 
sical papers by Velikhov [H and Chandrasekhar that 
a differentially rotating flow with a negative angular ve- 
locity gradient and a weak magnetic field is unstable to 
the magnetorotational instability. This instability has 
been analyzed in detail in the astrophysical context (see 
S 0) because it can be responsible for transport of 
the angular momentum in various objects ranging from 
accretion disks to galaxies. In accretion disks, this in- 
stability is also well studied by numerical simulations in 
both linear and non-linear regimes. Simulations of this 
instability in accretion disks (see, e.g., [1,0, Si) show that 
the generated turbulence can enhance substantially the 
angular momentum transport. 

Astrophysical applications of the magnetorotational 
instability rise great interest in trying to study this insta- 
bility in laboratory U, 12|. The experiments, how- 
ever, are complicated because very large rotation rates 
should be achieved. Recently, HoUerbach and Riidiger 
[ist argued that the rotation rate can be substantially 
decreased adding an azimuthal field. It is known since 
the paper by Tayler 1J| that an azimuthal field produce 
a strong destabilizing effect and, as a result of this ad- 
ditional destabilization, the critical Reynolds number in 
experiment can be reduced. 

On the other hand, the magnetorotational instabil- 
ity is not the only instability that operates in differen- 
tially rotating magnetized flows. For example, even a 
weak axial dependence of the angular velocity can re- 
sult in a double diffusive instability that is often called 
the Goldreich-Schubert-Fricke instability (e.g., [3], [3]). 
Note that many previous stability analyses have adopted 



the Boussinesq approximation, and have therefore ne- 
glected the effect of compressibility. This is allowed if 
the magnetic field strength is essentially subthermal, and 
the sound speed is much greater than the Alfven veloc- 
ity, Cs ^ CA but often this cannot be realized in real 
astrophysical conditions and in many numerical simula- 
tions. An attempt to consider the effect of compressibil- 
ity on the magnetorotational instability was undertaken 
by Blaes and Balbus [l3] in the context of astrophysical 
disks. The authors considered a very simplified case of 
the wavevector parallel to the rotation axis and a van- 
ishing radial magnetic field. As a result, the most inter- 
esting physics has been lost in this study since only the 
standard magnetorotational instability operates in this 
simple geometry. 

In this paper, we show that a new instability different 
from the standard magnetorotational instability may oc- 
cur in a compressible differentially rotating magnetized 
flow. This instability appears for any differential rotation 
and may occurs if the magnetic field has non-vanishing 
radial and azimuthal components. The instability can 
arise even in a sufficiently strong magnetic field that sup- 
presses the magnetorotational instability. Stability anal- 
ysis done in this paper will hopefully prove to be a useful 
guide understand various numerical simulations that ex- 
plore the nonlinear development of instabilities and their 
effects on the resulting turbulent state of rotating mag- 
netized flows. 



BASIC EQUATIONS AND DISPERSION 
RELATION 



We work in cylindrical coordinates (s, ip, z) with the 
unit vectors (e^, e^p, e^). The equations of compressible 
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MHD read 



il+{v-V)v = -—+g+-^{V X B)x B, (1) 



P + V-(pi;)=0, (2) 
p + v- Vp + 7pV • w = 0, (3) 

S - V X (u X B) + r?V X (V X B) = 0, (4) 

V-S = 0. (5) 

Our notation is as follows: p and v are the density and 
fluid velocity, respectively; p is the gas pressure; g is grav- 
ity that can be important in astrophysical applications; 
B is the magnetic field, r} is the magnetic difi^usivity, and 
7 is the adiabatic index. For the sake of simplicity, the 
flow is assumed to be isothermal. 

The basic state on which the stability analysis is per- 
formed is assumed to be quasi-stationary with the angu- 
lar velocity = ri(s) and B ^ 0. Generally, a quasi- 
stationary basic state cannot be achieved for any dif- 
ferentially rotating magnetic configuration, therefore we 
discuss in more detail when this assumption can be sat- 
isfied. We assume that gas is in hydrostatic equilibrium 
in the basic state, then 



— =D + —(V xB)x B , 
p Anp 



D = g + n'^s. (6) 



Generally, a geometry of the magnetic field can be rather 
complex, and our study is primcrily addressed such com- 
plex magnetic configurations where the radial and az- 
imuthal field components are presented. The presence 
of a radial magnetic field and differential rotation in the 
basic state can lead to the development of the azimuthal 
field. Nevertheless, the basic state can be considered in 
some cases as quasi-stationary despite the development 
of the toroidal field. 

For example, if the magnetic Reynolds number is large 
(77 is small), then one can obtain from Eq. (4) that the 
azimuthal field grows approximately linearly with time. 



B^{t) = B^{0)+sn'B,t, 



(7) 



where H' — dfl/ds, and -Bip(O) is the azimuthal field at 
t = 0. As long as the second term on the r.h.s. is small 
compared to the first one, and 



1 B^{0) 



(8) 



stretching of the azimuthal field does not aff'ect signifi- 
cantly the basic state; is the characteristic timcscale 
of generation of B^. As a result, the basic state can 



be treated as quasi-stationary during the time t <^ t^. If 
B^{0)/Bs 1, then steady-state can be maintained dur- 
ing a relatively long time before the generated azimuthal 
field begins to influence the basic state. We will show 
that the growth time of instability can be shorter than 
in many cases of interest. 

If the magnetic Reynolds number is moderate, then 
stretching of the azimuthal field from Bg by differential 
rotation can be compensated by ohmic dissipation, and 

the basic state can be quasi-stationary as well. Then, we 
have from Eq. (4) the following condition of steady-state 



or 



[V X (v X B)]^ = 77[V X (V X B)]^, 



A-\]B^ = --n'B,. 



(9) 



(10) 



The generated toroidal field is typicaly stronger than 
the radial field by a factor of the order of the magnetic 
Reynolds number. This simple model applies only in the 
case of moderate Reynolds number since the generation 
of a very strong toroidal field could lead to instabilities 
of the basic state caused, for example, by magnetic buoy- 
ancy or reconnection. Note that a quasi-stationary ba- 
sic state with non-vanishing radial and azimuthal field 
components can be achieved in other models as well. 
For example, if the angular velocity depends on both 
the s- and ^-coordinates, then changes in B^ caused 
by stretching from the radial and vertical field compo- 
nents due to radial and verical shear, respectively, can 
balance each other in such a way that will be steady- 
state. In fact, there is no principle difference for insta- 
bility which mechanism is responsible for maintaining a 
quasi-stationary basic configuration. The only impor- 
tant point for our model is the pesence of the magnetic 
field with non-vanishing radial and azimuthal compo- 
nents, but such magnetic configurations are rather com- 
mon in astrophysics (galactic and accretion disks, stellar 
radiative zones, oceans of accreting neutron stars, etc.) 

We consider the stability of axisymmetric short wave- 
length perturbations with the spacetime dependence oc 
exp{<Tt — ik ■ r). where k = {ks,0,kz) is the wavevec- 
tor, \k ■ r\ ^ 1. Small perturbations will be indicated 
by subscript 1, while unperturbed quantities will have 
no subscript. Then, to the lowest order in |A; • f|~^ the 
linearized MHD-equations read 



ikpi 



—(kx Bi) X B, 

A-Kp 



api — ip{k- vi) = 0, 



api — ijp{k ■ vi) =0, 



(11) 



(12) 



(13) 
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aBi = e^sn'Bu - i{B ■ k)vi + iB{k ■ i/i), (14) 

k -131=0. (15) 

We neglect ohmic dissipation in the induction equation 
because the inverse ohmic decay timescale is small for 
many cases of interest in astrophysics. 

Generally, the dispersion relation for Eqs.(ll)-(15) is 
rather complicated and, in this paper, we consider only 
a particular case when the wavevector of perturbations 
is perpendicular to B, k ■ B = Q. This case, being math- 
ematically much simpler, illustrates very well the main 
qualitative features of the new magnetic shear-driven in- 
stability. Besides, the standard magnetorotational insta- 
bility does not operate in this case because its growth 
rate is proportional to k ■ B. Therefore, the difference 
between instabilities is seen most clearly if fc • i? = 0. 

In the case k ■ B = 0, Eqs. (11)-(15) may be combined 
after some algebra into a fifth-order dispersion relation, 

a^+a'^iivl + nD + a^ujli^+annlijl+nnlijln = (16) 

where we denote 

nl = 2n{2n + sil') , u;', = k\cl + cl,} , fi = kl/k\ 
2 _ B"^ 2 _1P 3 _ k'^B^BsSn' 
VKp p Anp 

This equation describes five non-trivial modes that exist 
in a rotating magnetized flow if fc • = 0. 

In the non-magnetic case, B = 0, Eq. (16) yields 

a^ + {u;l + nl)a'+fiu;',nl = 0, (17) 

where Us = kcg is the frequency of sound waves. The 
solution is 

^2 = + ^l) ± + (is) 

Instability arises only if the well-known Rayleigh criterion 
is fulfilled, fi^ < 0. In this case, the inertial mode is 
unstable which corresponds to the upper sign. The sound 
mode that corresponds to the lower sign is always stable. 

To have an idea about the properties of dispersion 
equation (16), we can consider a particular case of flow 
with Vl (X s^^. Then, fig = and we have from Eq. (16) 

+ (7^2 ^ ^ 0. (19) 

The solutions of this equation are 

1 iVS 
ai=u + v, a2,3 = -^{u + v) ± —^{u - v), (20) 

where 



{u,v) = 




At least, one of the roots has a positive real part (insta- 
bility) ii u + V 0. The latter condition is equivalent 
^Bn 7^ ^ that is the criterion of instability in this simple 
case. It is clear from this simple example that the quan- 
tity LUBn plays a crucial role for stability of magnetized 
compressible flows. 

CRITERIA AND GROWTH RATE OF 
INSTABILITY 

The conditions under which Eq. (16) has unstable so- 
lutions can be obtained by making use of the Routh- 
Hurwitz theorem (see ^J). In the case of the dis- 
persion equation of a fifth order, the Routh-Hurwitz cri- 
teria are written, for example, in ^2d\. According to these 
criteria, Eq. (16) has unstable solutions if one of the fol- 
lowing inequalities is fulfilled 

pnlio%^<0, c^|o>0, (c^|n)'<0. (22) 

These inequalities yield the criterion of instability 

^ 0. (23) 

Apart from differential rotation, this criterion requires 
non- vanishing radial and azimuthal field components. 
The vertical component of B is unimportant for crite- 
rion (23), and the instability may occur even in a plane 
parallel magnetic field with components only in radius 
and azimuth. The direction of B and the sign of fl' are 
insignificant, and the instability may occur for both the 
inward and outward decreasing angular velocity. Note 
that this is in contrast with the magnetorotational insta- 
bility that can arise only if f2' < 0. Another important 
difference is that the magnetorotational instability is sup- 
pressed by a sufficiently strong field, whereas the insta- 
bility given by Eq. (23) can arise even in a very strong 
field. 

To calculate the growth rate in the general case it is 
convenient to introduce dimensionless quantities 

(we assume that fig > 0). Note that the parameters ^ 
and C do not depend on the wavevector. Then, Eq. (16) 
becomes 

This equation was solved numerically for different p, ^, 
and C by computing the eigenvalues of the matrix whose 
characteristic polynomial is given by Eq. (16) (see [2]| . 
for details). 

In Fig. 1, we plot the dependence of the real and imag- 
inary parts of r on a; for p = 0.3, ^ = 0.1 and = 0.1. 
The solid lines show the growth rate and frequency for 
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complex conjugate roots, and the dashed line for a real 
root. As mentioned, there should be no instability in the 
incompressible limit because all the considered perturba- 
tions are stable with respect to the standard magnetoro- 
tational instability. Our calculations, however, clearly 
indicate that some roots have a positive real part and, 
hence, there should exist a new shear-driven instability 
in the compressible flow with C 7^ 0. There are two pairs 
of unstable complex conjugate roots and one real stable 
root with negative Re T. In the considered domain of pa- 
rameters, Im r for complex roots is typically ~ 10 — 30 
times greater than Re F except the region of not very 
large ~ 10 — 50 where they are of the same order of 
magnitude (but still Im F > Re F). One pair of unstable 
roots has a very small growth rate ^ 10~^fle, but an- 
other one grows much faster. For these roots, the growth 
rate is ~ 0.5f2e and varies very slowly with the wave- 
length of perturbations. Note that calculations for other 
values of the parameters show that typically Re F ~ 0.5 
if C ~ C: but Re F becomes smaller ii (, ^ This is 
qualitatively clear because the case ^ C corresponds to 
the incompressible limit when the considered instability 
is substantially suppressed. 

In Fig. 2, we plot the dependence of Re F on a; for the 
case fj, = 0.3, ^ = 0.1 and C = —0.1. We do not plot Im F 
because this dependence does not differ much from what 
is shown in Fig. 1 . The change of sign alters qualitatively 
the behavior of roots. If C is negative then all oscillatory 
modes are stable (Re F < 0) but the real mode becomes 
unstable. This conclusion is completely consistent with 
our analytical consideration of Eq. (14). It is worth men- 
tioning that calculations for other /i, f , and C also indi- 
cate that this sort of behavior is rather general, and the 
non-oscillatory mode is typically unstable for negative ^ 
whereas the oscillatory modes are unstable for positive C. 
Like the previous case, the growth rate weakly depends 
on the wavelength except the region x"^ < 200 where this 
dependence is stronger. The characteristic value of the 
growth rate is larger for negative C and reaches « fte for 
> 200. Note that this is typical also for other values 
of fi and 1^ and that a non-oscillatory mode (negative C) 
grows faster than oscillatory modes (positive () for the 
same 

Fig. 3 illustrates the behavior of roots as functions of 
the parameter ( for fixed value oi x. It is seen that Re F 
vanishes for both oscillatory and non-oscillatory modes 
when C goes to zero. Since C cx ajsn, the instability 
occurs only if lobq ^ in complete agreement with the 
criterion (13). As usual, the real root is positive (instabil- 
ity) at C < whereas the oscillatory roots have positive 
real parts at C > 0. For the same |C|, the growth rate is 
larger for negative C.- 
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FIG. 1: The dependence of the real and imaginary parts of F 
on for fi = 0.3, ^ = 0.1, and = 0.1. Solids lines show the 
growth rate and frequency of complex roots, and the dashed 
line corresponds to the real root. 
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FIG. 2: The dependence of the growth rate on 3? for /i — 0.3, 
^ = 0.1, and C, = —0.1. The solid and dashed lines correspond 
to complex and real roots, respectively. 

DISCUSSION 

To summarize then, we have considered the instabil- 
ity caused by differential rotation of compressible mag- 
netized gas. To illustrate the main qualitative features 
of the instability associated to compressibility and shear, 
we analyzed a particular case of perturbations with the 
wavevector k perpendicular to the magnetic field B. 
In this case, the standard magnetorotational instability, 
well-studied in incompressible fluids (see, e.g., [l[, [1], 
[3]), does not occur because its growth rate is propor- 
tional to (k ■ B). Nevertheless, even perturbations with 
k-B = turn out to be instable if the necessary condition 
of the new instability, ^ oc BgB^Vi' ^ 0, is satisfied. 

In our stability analysis, we assume that the basic state 
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FIG. 3: The dependence of Re T on C for fi = 0.3, = 10 
and C = 1. Solid and dashed lines correspond to complex and 
real roots, respectively. 



is quasi-stationary. This assumption can be fulfilled in 
many cases of astrophysical interest despite the devel- 
opment of the azimuthal field from the radial one due 
to differential rotation. For instance, if the magnetic 
Reynolds number is large then the timescale of gener- 
ation of the toroidal field is ~ if B^{0) > Bg. In- 
stability can be considered in a quasi-stationary approx- 
imation if its growth time is shorter than t^. As it is 
seen from Eq. (21), the growth rate of instability in the 
case of a strong compressibility can be roughly estimated 
as LOsn- Then, the condition of quasi-stationarity reads 
^so > ^/t^p, or 



B^iO) 



(26) 



Since the l.h.s. of this equation is proportional to k but 
the r.h.s. does not depend on k, there always exists the 
range of k for which Eq. (26) can be satisfied and the 
basic state is quasi-stationary. 

The considered instability is related basically to shear 
and compressibility of a magnetized gas. In the incom- 



pressible limit that corresponds to Cg 
we have from Eq. (16) 



oo, 



(27) 



The properties of the considered instability are very 
much different from those of other instabilities which 
can occur in cylindrical magnetized flows. The neces- 
sary condition of the new instability (23) can be sat- 
isfied for both outward increasing and decreasing ri(s), 
whereas the magnetorotational instability occurs only if 
n{s) decreases with s. The found instability operates 
only if the basic magnetic configuration is relatively com- 
plex with non-vanishing radial and azimuthal field com- 
ponents while the standard magnetorotational instability 
can arise also if both these components are vanishing and 
only B, ^ 0. 

This new instability can be either oscillatory or non- 
oscillatory, depending on the sign of (, whereas the 
standard magnetorotational instability is always non- 
oscillatory. Typically, the considered instability is non- 
oscillatory if < and oscillatory if ^ > 0. One 
more important difference is associated with the depen- 
dence on the magnetic field strength. A sufhciently 
strong magnetic field, satisfying the inequality {k ■ B)^ > 
87rpsf2|f2'|(fc^//c^)|, completely suppresses the standard 
magnetorotational instability. On the contrary, the in- 
stability discovered in our study cannot be suppressed 
even in very strong magnetic fields as it is seen from the 
criterion (23). All this comparison allows us to claim 
that our analysis demonstrates the presence of the new 
instability in compressible cylindrical flow. 

The growth rate of the newly found instability can be 
rather large and reach ^ fig. Basically, the growth rate 
is larger for non-oscillatory modes which are unstable if 
^'bq. < 0. The growth rate depends on compressibil- 
ity, being smaller for a low compressibility. The incom- 
pressible limit (Boussinesq approximation) corresponds 
to Cg 3> CA, and the considered instability is inefficient 
in this limit because of a low growth rate. However, in 
the case of a strong field with ca ~ Cs when the Boussi- 
nesq approximation does not apply, the instability can be 
much more efficient than the magnetorotational instabil- 
ity. 
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and the instability disappears. This is a principle differ- 
ence to other well-known instabilities caused by differen- 
tial rotation such as the Rayleigh or magnetorotational 
instabilities. Note that an attempt to consider instabil- 
ity associated to compressibility of differentially rotating 
magnetized gas has been undertaken by Blaes & Balbus 
17| . These authors, however, analyzed only the unper- 
turbed configuration where the magnetic field has verti- 
cal or azimuthal components, but such configurations are 
stable in accordance to our criterion (23). 
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